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^ I The vicious random walker problem on a one dimensional lattice is considered. Many 

walkers take simultaneous steps on the lattice and the configurations in which two of 
them arrive at the same site are prohibited. It is known that the probability distribution 
of A^ walkers after M steps can be written in a determinant form. Using an integration 

^ \ technique borrowed from the theory of random matrices, we show that arbitrary fc-th 

I \ order correlation functions of the walkers can be expressed as quaternion determinants 

rH ■ whose elements are compactly expressed in terms of symmetric Hahn polynomials. 
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1 Introduction 

The vicious walker problem first introduced by Fisher [jl| and then developed by Forrester 
(see also ||^) recently attracts much attention in mathematical physics. Fascinating con- 
nections to other research fields, such as Young tableaux in combinatorics^, Kardar- 
Parisi-Zhang (KPZ) universality in the theory of growth process^ and the theory of 
random matrices^, have been revealed one after another. In the context of random ma- 
trix theory, the ensemble of vicious walkers in one dimension corresponds to a discretiza- 
tion of the Gaussian ensembles of random matrices. Therefore the theory of discretized 
random matrices is expected to shed light on all of the related problems. Indeed dis- 
cretized random matrices in the guise of discrete Coulomb gases are central to the work 
of Johansson m ^ |[ in these directions. 

Suppose that there are A^ walkers on a one dimensional lattice. In the lock step 
version of the model (as distinct from the random turns version considered in the recent 
work0), at each time step each walker moves to the left or right one lattice site with 
equal probability. Walkers are "vicious" so that two or more walkers are prohibited to 
arrive at the same site simultaneously. The j-th walker starts at the position x = 2j — 2 
and, after M steps, arrives at x = Xj. The walker configurations form nonintersecting 
paths in the x-t plane. An example is given in Figure |1|. Furthermore it has long been 
realized fiy] that there is a one-to-one correspondence between such nonintersecting paths 



and random rhombus tilings of a hexagon, or suitable truncation thereof, involving three 
types of rhombi. This is also illustrated in Figure |l|. 



The number of lock step paths is known to be expressed as a binomial determinant |11 
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The binomial determinant can be further rewritten as a product formula 
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(1.2) 
where Xi < X2 < ■ ■ ■ < X^- Now we introduce new variables 



Xj 



Xj-N +1 



L = M + N-1 (1.3) 

in order to obtain a compact and symmetric expression 

N N 

Pi{xi,X2, ■ ■ ■ , xn) = Cmn n y'^i^j) n I ^i ~ ^' I • (1-4) 



Here 



N 



Cmn = l[(M + N - j)] (1-5) 





Figure 1: Nonintersecting paths representing a vicious walker configuration with A^ = 3 
walkers and M = 4 steps. In the top diagram these paths are drawn in the x-t plane. In 
the bottom diagram superimposed on the paths is the equivalent rhombi tiling, involving 
left sloping rhombi (step to the left), right sloping rhombi (step to the right) and vertical 
rhombi. 
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In some applications of the vicious walker problem one imposes the additional con- 
straint that each walker returns to the initial position after 2M steps. One example is the 
rhombus tiling problem of Figure [l| with the region to be tiled extended to be symmetrical 
about t = M and thus made into a hexagon|^. In such cases the number of paths is given 
by the square of Pi(xi, X2, ■ ■ ■ , xat): 



P2{Xi,X2, 



Xn 
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As will be revised below, the function w{x) is a special case of the Hahn weight function 
from the theory of discrete orthogonal polynomials. The probability density ( p..7|) with 
the Hahn weight is intimately related to Hahn polynomials and so has been termed the 
Hahn ensemble 0, ||. 

We are interested in the number of paths under the condition that k walkers take fixed 
positions Xi, X2,- ■ ■ ,Xk after M steps {k < N). This number is given by the correlation 
functions 



r(/3) 



II lxi,a;2. 
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in both the cases P = I and P = 2. 
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In the case with returning walkers {f3 = 2), the evaluation of the correlation func- 
tions is relatively easy. Introducing nionic orthogonal polynomials Cj (x) = x^ + ■ ■ ■ with 
orthogonality relations 

oo 

Y, wix)C,{x)Ci{x) = 5jihj, (1.9) 

a;=— oo 

we can readily obtain 

N-l 

ll^\xi,X2r--,Xk) = Cl^N n hj det[K{xj, xi)]j^i=i^2,-,k (1-10) 

3=0 

with 
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K{x,y) = ^w{x)w{y)J2 —Cj{x)Cj{y) 

j=o '^i 

1 CM{x)CM-i{y)-CM-i{x)CM{y) 



w{x)w{y)- . (1.11) 
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We call Cj{x) symmetric Hahn polynomials and define them in terms of the Hahn poly- 
nomials in §2. 

The main purpose of this paper is to derive an analogous formula for the correlation 
functions with no returning constraint {(3 = 1). For that purpose, we again make use of 
the above symmetric Hahn polynomials to rewrite the correlation functions in the form 
of a determinant although now with quaternion elements. The Hahn polynomials and 
their symmetrization are introduced in §2. In §3, we introduce quaternion determinant 
formulas for the correlation functions I^ (xi, ■ ■ ■ ,Xk)- In §4, the continuous limit to the 
Gaussian ensembles is discussed. 



2 Symmetric Hahn Polynomials 

The Hahn polynomials Qn{x) are orthogonal polynomials with respect to the weight 
function 1 1^ 

{a + x)\{L + b-x)\ 

"-^"^ = alxl U{L-x)\ ^'-'^ 

on a discrete measure (x is restricted to be an integer). When the coefficient of the highest 
order term is fixed as 

(a + l)n[-L)n 



we have the orthogonality relation 

L 
Y^ WH{x)Qm{x)Qn{x) = HnSmn, (2.3) 
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where 



" L\{2n + a + b+l){-L)n{a + l)n' ^'' 



Here {a)n = T{a + n)/r{a). Moreover there is a dual orthogonality relation (completeness 
relation) 

E ^Qnix)QM = -%T, (2.5) 

n=0 " ^H\-^) 

which yields the orthogonality relation of the dual Hahn polynomials. 

The weight ( |1.6| ) is the special case a = b = —L — 1 of ( p.l| ) with shift of coordinate 
X I— > a; + L/2. The monic orthogonal polynomials with respect to (|1.6|) , which we have 
denoted C„(x) and termed the symmetric Hahn polynomials, are therefore given in terms 
of the Hahn polynomials Qn{x) by 



Cn{x) — - ; I , , tn Qn ( ^^ + — 

{n + a + b + l)n V 2 



(2.6) 



a=b=-L-l 



and satisfy the orthogonality relation 

E w{x)C,{x)Q{x) = S,ih, (2.7) 

x=—oo 

with 

_ j!(2L-2j + l)!(2L-2j)! 

^■" (2L-j + l)![(L-j)!]^ • ^ ^ ^ 

The corresponding dual orthogonality relation reads 

±l.C„(.)C„(y) ^ ^^. ,2.9) 



Since the weight function w{x) defined in (pTq ) is symmetric (w(— a;) = w{x)), Cn{x) are 
even or odd polynomials corresponding to the parity oi n (C„(— x) = (— l)"C„(a;)). 



In Ref. |15|, recurrence relations for the Hahn polynomials are presented. In the sym- 



metric limit a = b = —L — 1, we find 

Cn+l{x) = xCn{x) - UJnCn~l{x) (2.10) 

and 

-(^+^) [Cn{x) - Cn{x - 1)] = a„C„+i(x) + ^nCn{x) + 7nC„_l(x), (2.11) 



where 



and 



n(2L-n + 2)(L-n + l)^ 
"^ 4(2L-2n + 3)(2L-2n + l)' ^ ' ' 

a„ = -n, /5„ = --n(2L - n + 1) (2.13) 

n(2L-n + l)(2L-n + 2)(L-r2 + l)^ 
^"~ 4(2L-2n + 3)(2L-2n + l) ' ^ ' 



3 Correlation Functions 



In the theory of random real symmetric matrices, it is known that the correlation functions 



among the eigenvalues are expressed as quaternion determinants [0, [T^, |T8[. Although in 
most cases the integration method has been applied to continuous eigenvalue distributions, 
several authors have made attempts to evaluate the correlation functions on discrete 



measures|T^, [2^, ^ ^. Here we can employ a similar procedure to rewrite the multiple 
sum (|l.(j|) in a quaternion determinant form. 



3.1 Quaternion Determinant 

As the first step let us introduce the quaternion determinant. A quaternion is defined as 
a linear combination of four basic units {1, ei, 62, 63}: 



g = go + q-e = go + giCi + ^262 + ^363, 



(3.1) 



where qo,qi, q2 and ^3 are real or complex numbers. The first part gi is called the scalar 
part of q. The multiplication laws of the four basic units are given by 



1-1 = 1, 1 ■ e,- 



e^-l 



^j' 



616263 



J = 1,2,3, 



-1. 



(3.2) 
(3.3) 



Note that the quaternion multiplication is associative but in general not commutative. 
The dual g of a quaternion q is defined as 



g = go - q • e. 



(3.4) 



A matrix Q with quaternion elements qji has a dual matrix Q = [qij]. The quaternion 
units are represented as 2 x 2 matrices 
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(3.5) 



so that any 2x2 matrix with complex elements can be identified as a 2 x 2 representation 
of a quaternion. 

We define a quaternion determinant Tdet of a self dual Q ( i.e., Q = Q ) as 



Tdet Q = E(-l) Ui^atqbc ■ ■ ■ qdaV 
p 1 



(3.6) 



Here P denotes any permutation of the indices (1, 2, ■ ■ • , A^) consisting of / exclusive cycles 
of the form {a^h^c^---^>-d^>-a) and (—1)^"' is the parity of P. The subscript 
has a meaning that we take the scalar part of the product over each cycle. 



3.2 Skew Orthogonal Polynomials 



Next we need to define skew orthogonal polynomials in order to give the quaternion 
determinant expressions for the correlation functions. The Schmidt's orthogonalization 
procedure enables us to construct monic polynomials -R„(x) of degree n which satisfy 



{R2m{x), R2n{y)) = 0, {R2m+lix), i?2n+l(2/)) = 0, 



(3.7) 



where 



{f{x),9iy)) = ^ E 

y=-oo 



y-1 oo 

5=— oo x=y+l 



\/^)9iy)- (3.8) 



When the weight function w{x) is given by eq. ( |1.6| ), the skew orthogonal polynomials 
Rn{x) are compactly expressed in terms of the symmetric Hahn polynomials. To show 
this we utihze a method analogous to that employed by Nagao and Wadatij^ for the 
continuous analogue of ( |3.8| ). Let us first note the identity 



{f{x),g{y + l)-g{y)) 



-{f{x)^ r^= 9{y)) 



w{y) 
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Y. \MaO/(^) yu!(x)g{x + 1) + \Jw{x - l)g{x) 



and make substitutions 



Noting 
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Ci{x) + \Jw{x) 
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2-'. 
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+ A 
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+ ^Mx)[{L-l)uJi + -ii]Ci^i{x), 
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(3.10) 
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we obtain 
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+ o Z^ w{x)Rn{x) j^^ o/ , o^/o^ ^^^^ Ci.i[x). 



Putting the expansion 



(2L-2/ + 3)(2L-2/ + l) 



R2m{.x) = Ya2m2jC2j{,x), 

j=0 

rn 

R2m+l{x) = 2^<y2m+l 2j+lC2j+l{x) 

j=0 

{ajj = 1) into ( ^.lij| ) and using eqs. ( |2.7] ) and (|3.7D, we arrive at 

R2m{x) = C2m{x), 
R2m+l{x) = C2m+l{x) "TT C2m-l{x) 



L-2m + lh 



2m-l 



and 



-m = -{L -2m)h2m, 



where hn is given by eq.(^ 



(3.13) 



(3.14) 



(3.15) 
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3.3 The Case N Even 

Let us first consider the case A^ is even (A^ = 2z/). Independent of the functional form of 
w{x), the number of paths has a quaternion determinant expression [O 



u-l 



Pl{Xi,X2, ■■■,Xn) = CMN2''iY[ rn)Tdet[f{Xj, Xi)]j^i=i^2,-,2u, 

n=0 



where the quaternion elements are represented as 

f{x,y) 

with 

AN/2)-l 



S{x,y) I{x,y) 
D{x,y) S{y,x) 



S{x,y) = yw{y) Y. —[^2n{x)R2n+l{y) - ^2n+l{x)R2n{y)], 
^ — n I n. 



n=0 



.(N/2)-l 



D{x,y) = Jw{x)w{y) J2 — [^2n(a;)i?2n+i(l/) - ^2n+i(a;)i?2n(l/)], 



n=0 ^iT- 



I{x,y) 



(7V/2)-l 



n=0 ^" 



(3.17) 



(3.18) 



J2 —[^2nix)^2n+liy) ~ $2n+l (a;)$2n(l/)] - e{x - y) (3.19) 



and 

OO 

^n{x)= E e{x-y)^w{y)R^{y). (3.20) 

y=-oo 

Here e(x) is defined as 

r 1/2, a; > 0, 

e{x) = \ 0, a; = 0, (3.21) 

[ -1/2, x<0. 

Using the skew orthogonality relation ( p.7| ), we can prove the following. Let us define 

Qn(a;i,X2,---,a:„) = Tdet[/(xj, X;)]j,i=i,2,...,n- (3.22) 

Then 

OO 

E Qn{Xi,X2r-- ,Xn) = {N -n+l)Qn-l{Xi,X2,---,Xn-l). (3.23) 

Recursive use of (|3.23|) on (|3.22|) leads to 



jf ^ (xi, X2, • • • , Xfc) = CMAr2''( n r„)Tdet [/(xj, Xi)]j,i=i,2,...,fe. (3.24) 

n=0 

Thus the correlation functions are expressed in a quaternion determinant form. 

In the continuous case, it is known [^ that for w{x) a classical weight the summations 
( p.l9| ) can be simplified. Similar simplifications can be undertaken here. First we note 
the function ^n{x) can be rewritten as 

1 OO OO 2' — 1 

^n{x) = - ^ ^ Yl l^yz^xz' -SxzSyz']yw{y)Rn{y). (3.25) 

■^ y = -QO z' = -oo z = -oo 



We put the dual orthogonality relation ( |2.9| ) and an expansion 

m 
C2m{x) = / ,/^2m 2jR2j{x), 



j=0 
m 

C2m+l{x) = / ,/^2m+l 2j+lR2j+l{x) (3.26) 

i=0 



in the above equation to find 



^2n{x) = \lw{x) Yl -T^ /^2p+l 2n+l ?"„, 

p=n l>'2p+l 

■^"'^C2p{x^ 



'^2n+l(a;) = -Jw{x) J2 —r (32p 2n rn- (3.27) 



p=n '^2p 



Here [x] is the largest integer not exceeding x. Substituting (|3.27| ) into (|3.19| ) and using 
( p7[5| ) yields 

I ^-^ 1 

S{x,y) = \Jw{x)w{y)Y T-Cn{x)Cn{y) 

n=0 '^" 



.[(L-l)/2](7V/2)-l 



+ ^W{x)w{y) Y. H -j;^ /^2p+l 2n+1^2n+l(y) 

p=N/2 n=0 '^2p+l 

I ^-^ 1 

= ^w{x)w{y) Y^ —Cn{x)Cn{y) 

n=0 '^" 

+ ^w{x)w{y)CM-i{y) E P2p+iN-i^^^^. (3.28) 

p=N/2 '^Sp+l 



Making use of 



^ / \ I TT v^ o C2p+i{x) I -—Cn~i[x) 

^N-2[X) = ^Jw[x) 2^ (32p+l N-l—r r(Ar/2)-l + \^ w[x)—- r(Ar/2)-l (3.29) 

p=N/2 "'2p+l rLN-1 

and i?7v-2(a^) = C]^^2{.x)i we can easily find 

Af-2 1 

S{x,y) = ^/w{x)w{y) J2 —Cn{x)Cn{y) 



71=0 ""n 



1 / °° I 

+ ^w{y)CM.M E e{x-z)^w{z)CN_2{z). (3.30) 

^{N/2)-l ^=_oo 

Note that the first term on the right hand side is precisely ( |1 . 1 Ij ) with N ^-^ N — 1. This 
gives a compact expression of S{x, y) in terms of the symmetric Hahn polynomials. 

It should be also noted that, using eq. ( p.27| ), we can derive 

[(i-l)/2] ^ 
e{y-x)= E —['^2nix)^2n+l{y) - ^2n+lix)<^2niy)l (3.31) 



which yields 



n=0 ''"n 



[(i-l)/2] ^ 



I{x,y) = E -['^2n(a:)$2n+i(y) - $2n+i(a;)$2n(y)]. (3.32) 



n=N/2 ^" 

To conclude this subsection, let us remark on the "partition function" 
4') = CA/7v2^(nr„) 

n=0 

^^/g)-^ (2n)!(2L-4n + l)!(2L-4n)! 
J-io 2(2L-2n + l)![(L-2n)!p ' ^ ' ^ 

It can be readily rewritten as 

(i)_^^ 1 ^^/^''^ {2n)\{2M + An + 3)\{2M + An + 2)\{2M + 2n + 1)\ 

° ~ /j (2M + 2/ + 1)! J-j, 2[(M + 2n + l)!]2 ' 

(3.34) 
so that 






Ar^;v+2_ N\{2M + N+l)\ _ ^ N + t+j + 1 



j(i) 
-'o 



{M + N + 1)\{M + N)\ ^^^^^^N + t+j-1 ^^-^^^ 
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Therefore, noting 



r(l) 



1 (2M + 2)! 



N^2 



2[(M + 1)!] 



„ i+j+l 



we obtain 



r(l) 



n 



i<«<i<Af 



A^ + ^ + J - 1 



(3.36) 



(3.37) 



This expression is aheady known. It was conjectured by Essani and Guttmann [|ri|] and 



subsequently proved in Ref. |12| using the correspondence between the paths of vicious 
walkers and Young tableaux. In the context of Young tableaux, this enumeration problem 
is reduced to a special case of the result known as Bender-Knuth conjecture. 



3.4 The Case N Odd 

In order to deal with the case N odd (A^ = 2i/+ 1), we introduce the polynomials R'^^{x) 



K^\x) = Rn{x) 



SN-1 



-Rn^i(x), n = 0,l,2,---,iV-2 



where 



^n= Yl \Jw{x)Rn{x). 



(3.38) 



(3.39) 



Then we know from ||2^ that we can express the number of paths in a quaternion deter- 
minant form 



v-l 



Pi(xi,X2, ■ ■ ■ ,x^) = Cmn sn^i 2'^(n r„,)Tdet[/°<^'^(a;,-,a;0],v=i,2,.,2.4-i, (3-40) 

n=0 

where the 2x2 representations of the quaternion elements are given by 



5°<^^(a;,2/) P'^'^{x,y) 
D°'^'^{x,y) S'"^'^{y,x) 



(3.41) 



with 



.{iV-3)/2 



S'"'\x,y) = ^) Y: -[<'(^Xti(2/)-<+i(^X'( 

+ Jw{y)RN-i{y), 

sn-i 

^ (JV~3)/2 . 

D^'^'ix^y) = ^w{x)w{y) ^ -[i?^ (xXt,(i/) - <t,(x)i?°f (y)]. 



n=0 ^ri- 



(7V-3)/2 



1 



n=0 ^" 



SN-1 



[^N-i{,x) - ^N-iiy)] - e(x - y) 



(3.42) 



11 



and 

OO 

<\x)= Y: 6(a:-y)vMy)<'(2/)- (3.43) 

y=-oo 

The skew orthogonality relation among Rn{x) enables us to carry out the multiple inte- 
gration one by one as before and arrive at 

u-l 

li'\xi,X2,---,Xk) = CMN SN-1 2^(n^n)Tdet[r'^'^(x,,xO],7=i,2,.,fc. (3.44) 

n=0 

This is the quaternion determinant formula for the correlation functions in the case N 
odd. 



Inserting (|3.38| ) into (p.42| ), we can rewrite S°'^'^{x,y) as 



sn-i 



+ ^-^ E —[<^2n+l{x)RN-l{y)-^N-l{x)R2n+l{y)]. (3.45) 
SN-1 n=0 ^n 



If we note from eqs. (|3.19|) and (|3.30|) that 



. {N-3)/2 

S2n 



w{y) J2 -^R2n+i{y) = 2\imS{x,y)\^^^^^ 



n=0 ^n 



^^^N^^C^^M, (3.46) 

r{N-2)/2 



this can be further simplified as 



S^^^x.y) = S{x,y)\^^^_, + ^^CM-i{y) 

sn~i 

SN-3 



+ yMy) ^^"^ [ipN~2{.x)CN-i{y)-^N-i{x)CN-2{y)], (3.47) 

SN-1 r(N-3)/2 



where 



¥'n(x)= E eix-y)^wiy)Cniy). (3.48) 



y=-oo 



We remark that s„ as defined by ( |3.39| ) can be written in a closed form. For this 



purpose we require the generating functions for the dual Hahn polynomials in Ref.|T3 
One of them can be written in terms of the Hahn polynomials as 



(1 - t)^-" 2F,{-n, -n-b;a + l;t) = J2 ^-^Qn{x)t^ 



(3.49) 



x=0 



Putting t = — 1 and a = b = —L — 1, we deduce 

. „n!(L-(n/2))!(2L-2n + l)! 



n even, 



(n/2)!(2L-n + l)![(L-n)!]2' ' (3.50) 

0, n odd. 
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As an application, since the "partition function" Iq can be rewritten as 
4'^ = Cmtv s^-i 2^(n r„) 

n=0 

^"^-^ ■ H (2M + 2/ + 3)! J-j, 2[(M + 2n + 2)!]2 ' 

(3.51) 



we can utilize the summation formula ( |3.50|) to derive 
r(i) 



I^ 



N^N+2_ N\{2M + N + 1)\ _ ^ N + t+j + 1 



I'o 



r(i) 
which leads to 



N^N 



(M + iV + l)!(M + iV)! ^^^ii^^^iV + ^+j-1' ^- ' 



4^^= n ,,,_i • (3-53) 

l<i<j<M '"""J -^ 



This is again a known result [|T2|, related to the Bender-Knuth conjecture in the theory 
of Young tableaux. 



4 Continuous Limit 

In this last section we discuss the continuous limit to the Gaussian ensembles. Some 
aspects of this limit for the /3 = 2 Hahn ensemble ( |1.7|) have been discussed in ^. Let us 
first define 

nL (D.y fy/L 



«^(^) = — I ^ 1 ^ I ^^^ 1 • (4-1) 



Then the Stirling's asymptotic formula for the Gamma function yields 

w^^\x) = lim w{x) = e""^". (4.2) 

L— ►oo 

Hence we can say that the vicious walker problem is reduced to the Gaussian ensembles 
of random matrices in this scaling limit, with ( |1.4| ) identical to the probability density for 
the Gaussian orthogonal ensemble. 

We further define monic polynomials 

CM=[^^ a(^x^ (4.3) 

and 

CPix) = lim Cn{x). (4.4) 

Then, from eg. ( ^.71 ), it can be readily seen that 
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In the asymptotic limit L ^ oo, the infinite sum in the above equation becomes a Riemann 
integral and results in 

J2^e-^'cf\x)C}''\x)dx = V^^^6,i, (4.6) 

which means 

CL'^Kx) = ^H4x), (4.7) 

where Hn{x) is the Hermite polynomial. The corresponding skew orthogonal polynomials 
are introduced as 

Rnix)=(-^) Rn(^x\ (4.8) 

and 

R^^\x) = lim Rn{x). (4.9) 

L— >oo 

Then we can similarly find 

{RS!(x),R^i\y))G = 0, (42+i(a;),4nli(2/))G = 0, (4.10) 

where 

{fix),giy))G = - dy dx e'^' /' f (x) - dx e-^'^'fix) e'y'/'giy). (4.11) 

Z J —oo L-^— oo Jy 

The scaling limit of eq. (p.l5|) 



R2m\X) — (~^2m\X), 

RS!+iix) = Ci%,{x) - mCS-,{x) (4.12) 

reproduces the known compact expression of Rl^\x) (see e.g. ||18|, ^, |4|)- 



The Gaussian ensembles have been thoroughly studied by many authors and asymp- 
totic properties in the limit N ^ oo are well known. In the vicious walker problem, 
Gaussian results are valid only when we first take the limit M — > oo in the above way 
and then analyze the large A^ asymptotic behavior. Other scaling limits with respect to 
M and N are not so well understood and should be further investigated in future works. 
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